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Quantification of red soil macropores affected by slope erosion and 
sediment using computed tomography 


Soil structure mechanism of collapsing gullies: Soil macropores, topography, 
water erosion and mass erosion 

Abstract: 

Purpose Soil structures are the main course of the formation and development of 
collapsing gullies, which are the most severe type of erosion in south China. 
However, few studies have focused on the relationship between soil macropores, 
soil erosion, and local topography. This study aimed to quantify and compare soil 
properties and macropore characteristics in the collapsing gully region, and 
explore their influences on the formation and development of the associated 
erosion. 

Materials and methods Soil core columns at different positions of a typical 
collapsing gully were excavated, and then scanned to analyze soil macropores. 
Moreover, soil properties and saturated hydraulic conductivity were investigated 
in the laboratory and in the field, respectively. 

Results and discussion The results indicated that the sand content increased from 
the ridge to the slope and the valley, while silt and clay contents decreased for the 
same catena. The mean weight diameter of aggregates was largest at the ridge and 
lowest at the valley. The infiltrate rates were highest at the valley and lowest at 


the slope. The valley had the greatest macroporosity (1.09+0.33%), and the 
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highest number (5919+703), volumes (24.7 3- 7.5 cm?) and surface (10.4 2.6 m?) 


of macropores, as well as the highest conectivity (42.3), while the slope had the 
smallest macroporosity (0.15+0.14%), and the smallest number (1189+747), 
volumes (3.4+3.2 cm?) and surface (1.71.4 m?) of macropores. The mean pore 
volume of macropores larger than 1 mm? was largest at the ridge (16.8+7.4), and 
smallest at the slope (10.6+2.9). The number of macropores and their 
macroporosity mainly decreased with increasing depth, but were influenced by 
the soil macrofauna as well as the erosion and sediment processes. Macropores 
were mainly vertical, which is affected by the roots of plants and is conducive to 
the vertical infiltration of water. But, there wer many horizontal macropores at 
the valley because of the sediment process. The equivalent pore diameter of 
macropores was mainly smaller than 2 mm (accouting for more than 76.3%), and 
the macropores larger than 5 mm were less than 196. 

Conclusions The macropore characteristics at different sites of the collapsing 
gullies affected the soil water infiltration and hydraulic conductivity, and further 
affected the processes of water erosion and mass ersion. The highest 
macroposities at the valley would result strong subsurface flow erosion and the 
loss of the base of collapsing wall. Macropores at the ridge would increase rain 
infiltration and promote soil collapsing. Few macropores and low infiltration 
abilities at the slope would strengthen the overland flow erosion. Thus, macropore 
characteristics had significant effects on both the formation and development of 


collapsing gullies. 
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collapsing gully 


1. Introduction 


Collapsing gullies, also called “Benggang” by locals, are the most 
extensively eroded areas of the seven provinces of South China. Both mass 
erosion and water erosion occur synergistically on the weathering crust of granite, 
at a depth of dozens to hundreds of meters, in tropical and subtropical regions. 
This area contains more than 2.4x10? collapsing gullies, covering an area of 
1.1x10° km? (Figure 1a), 45% of which are in Guangdong Province, which 
contains 74% of their total area (Liang et al., 2009). More than 6.7x10/ t yr! of 
soil is lost because of these collapsing gullies, which damages farmland, houses, 
roads, bridges, reservoirs, and ponds, and causes economic losses of 3.2x10? 
USD, affecting 9.2x10° residents (Liang et al, 2009). The occurring and 
developing mechanisms of collapsing gullies have rose many efforts of 
researchers, and soil macropore and soil water movement were thought to be 
important factors (Tao et al., 2020). 

Many factors influence both the spatial distribution and movement of soil 
water and other soil hydrological processes (Ghasemizade and Schirmer, 2013). 
The soil properties that particularly affect the soil hydrological process include 


soil structures, particle size distributions, porosities, soil layer characteristics, 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


90 


91 


p» 


93 


94 


95 


96 


97 


98 


99 


oxidation-reduction characteristics, organic matter contents, soil water repellency, 
and the soil water content itself (Zhu et al., 2012). Among these factors, the 
influence of macropores on soil water movement is of key importance (Lin et al., 
2010). Macropores allow rainfall to infiltrate into deeper soil layers, while 
impervious layers without macropores, such as caliche, prevent the continuous 
infiltration of water, thus leading to a higher soil water content above the caliche 
(Zhang et al., 2012; Li et al., 2013; Hu et al., 2018). The macropore structures of 
soil in different slope positions differ because of the different characteristics of 
soil aggregates, gravel contents, and plant roots. This results in differences in soil 
water conductivity and preferential flow characteristics (Holden, 2009; Zhang et 
al., 2016; Mei et al., 2018; Hu et al., 2020a). However, the preferential flow forms 
different macropore structures because of erosion (Nguyen et al., 2019). 

The tropical and subtropical regions of South China are subject to a lot of 
storm, but the alternations between dry and wet are obvious. Water is the limiting 
factor of the ecosystem service function during the dry season and also the driving 
factor of collapses, landslides, debris flows, and other disasters during the rainy 
season (Liang et al., 2014; Tang et al., 2014). The granite weathering crust in the 
study area can be divided into a topsoil layer, red soil layer, sand layer, debris 
layer, and spherical weathering layer from the top to the bottom (Tao et al., 2017). 
It has been reported that the infiltration capacity of the red soil layer, sand layer, 
and debris layer decrease sequentially from surface to the deeper soil layers 
(Duan et al., 2018). The water holding capacity of the topsoil layer and the red 
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soil layer is stronger than those of the sand layer and the debris layer. The soil 
water content gradually increases from the top to the bottom within the topsoil 
layer. However, the soil water content of the red soil layer to the debris layer 
gradually decreases. The farther away from the collapsing wall the site is, the 
higher its soil water content will be (Yusong ef al., 2015; Ni, 2016). The 
preferential flow varies by slope positions and soil layers (Tao et al., 2017). There 
are more preferential flows at the red soil layer than at the sandy soil layer and 
the clastic layer (Zhao, 2016), which are affected by soil properties, precipitation 
characteristics, and the initial water content (Duan et al., 2016). At present, 
studies of soil hydrology of collapsing gullies focus mostly on the characteristics 
and processes of hydrology, or on the characteristics of the physical and chemical 
properties of the soil. However, little research focuses on the comprehensive 
effects of geomorphology, vegetation, and soil properties, or their influence 
mechanism on soil hydrological process. Therefore, the influence of the soil 
macropore structure on soil hydrological processes, the mechanism of spatial 
distribution and transport of soil water, and their responses to future climate 
change are not fully clear. 

The loose and deep granite crust forms the material basis of collapsing 
gullies, and both rainstorms and soil water accumulation are the main drivers of 
the development of collapsing gullies. Soil structural properties (e.g., soil particle 
compositions, bulk densities, and porosities) significantly vary in different layers 
of the deep granite crust (Liao et al., 2018). This variation results in greatly 
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differing soil mechanical characteristics in response to wetting (Xia et al., 2018). 
The increase of the soil water content after rainfall and infiltration will lead to an 
increase of upper soil weight and a decrease of soil shear resistance, which in turn 
results in the collapse of the granite soil wall (Zhang and Zhong, 1990). At the 
same time, the decrease of the cohesion and internal friction angle of the granite 
soil, as well as the increase of the pressure of the soil fissure flow when subjected 
to water are the main causes of the granite soil wall collapse (Li, 1992). 
Preferential flow along soil macropores is considered to be an important reason 
for the spatial-temporal variation of the soil water content and soil fissure water 
movement (Hencher and Beven, 2010). 

Therefore, it is of great practical importance to study the soil macropores 
of granite red soil, which will aid investigations of the hydrpedological process 
toward a better understanding of the occurrence and development law of 
collapsing gullies. This study also helps to effectively prevent, control, and 


restore erosion-related damage in these regions. 


2. Materials and Methods 


2.1. Study sites 


The soil cores were sampled at a collapsing gully at Huacheng Town 
(24°04'N, 115°38'E), north of Wuhua County, Guangdong Province, China 


(Figure 1). Collapsing gullies are often close to residential and agricultural 
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regions, and the often devastating soil erosion and mud sand flow events (Zhang 
and Liu, 2014) increase the risk for houses, roads, and farmlands (Figure 1b). 
Wuhua County is a representative subtropical monsoon region with a mean 
annual air temperature of 21.3 °C, and a mean annual precipitation of 1507.2 mm 
(Zhang et al., 2020). The soil parent material is biotite granite and the weathering 
crusts are very deep, often exceeding 200 m. The soil is barren sandy loam with 
only 1.71 g kg! of soil organic carbon (Zhang et al., 2020). The vegetation is 
secondary with a coverage varying from 20% to 80% under different topographic 
and erosion conditions (Figure 1b). Pinus massoniana is the dominant 
constructive species, while Dicranopteris dichotoma, Rhodomyrtus tomentosa, 
Miscanthus floridulus, and Baechea frutescens form the main undergrowth 
(Zhang et al., 2020). Three sites along a slope near a Benggang, the ridge site, 
slope site, and the valley site (Figure 1b), were chosen and a total of nine 
undisturbed soil cores were collected in three duplicates per site. The ridge and 
slope sites were covered by sparse P. massoniana with rare undergrowth 
vegetation, whereas, the valley site was mainly covered by Phyllostachys edulis 
with dense D. dichotoma and Nephrolepis auriculata as undergrowth. 
«Figure 1> 


2.2 Analysis of basic soil properties 


For each site, three replicated soil samples were collected and soil 
properties were analyzed. Soils were air-dried at room temperature before sieving 


to 2 mm to analyze the soil particle size distribution with the pipette method (Gee 
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and Bauder, 1986). The aggregate size distribution was measured by separating 
different aggregates via wet sieving following a previously published method 
(Dong, 1997). First, field-moist soil samples were carefully passed through an 8- 
mm sieve (without breaking the soil structure) and then air-dried. Then, about 50 
g of air-dried soil for each replicate was wet-sieved through a series of four sieves, 
which separated the samples into five different aggregate size fractions: >2000 
um large macroaggregates, 1000-2000 um middle macroaggregates, 500-1000 
um small macroaggregates, 250-500 um tiny macroaggregates, and 106-250 um 
microaggregates. The four sieves were stacked so that the largest was at the top 
and the smallest at the bottom, and the soil sample was put into the top 2-mm 
sieve. Before wet-sieving, the soil samples were immersed in water on the top 
sieve for 5 min to break down aggregates via slaking pressure. Sieving was 
controlled by an Aggregate Analyzer (QT-WSI021, Qudaotech, Shanghai, China) 
with an up and down movement frequency of 30 times per min for a sieving 
period of 2 min. After sieving, each sieve was backwashed and the fraction that 
remained on the top of the sieve was collected, oven-dried at 105 °C, and weighed. 
The mean weight diameter (MWD) was computed as the sum of the oven-dried 
mass fraction of the soil remaining in each sieve after sieving, multiplied by the 


mean diameter of adjacent meshes (Le Bissonnais, 1996). 
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2.3 Field-saturated soil hydraulic conductivity measurements by DualHead 
Infiltrometers 


The field-saturated soil hydraulic conductivity (Ks) was measured using a 
DualHead Infiltrometer (Decagon Devices, Pullman, WA, USA) following an 
optimized method of Zhang et al. (2019). Before data collection, grass and 
vegetative litters were carefully removed and the surface soil temperature and 
water temperature were measured using ECH2O STE sensors (Meter Group, 
Pullman, WA, USA). An insertion ring with a depth of 5 cm and a radius of 7.5 
cm was gently driven into the soil to ensure good contact with the soil and 
minimize disturbance. It was ensured that the ring was level in all orthogonal 
directions. Ky was measured using a single cycle with 35 min holding time at 
each pressure head of 10 and 15 cm after 15 min of soaking. Only measurements 
without leakage and with coefficients of variation of the pressure head and water 
level below 15% were assumed to be reliable and used for further analyses. All 
Ky; values were viscosity-corrected to a standard temperature of 25 °C before 
comparison among different sites. This correction prevented the influence of 
changes in water effluent viscosity at different temperatures according to a 
previously described method (Zhang et al., 2019). 


2.4 Soil core samples 


Polyvinyl chloride (PVC) cylinders with 5-mm wall thickness, 10-cm 
diameter, and 50-cm length were used to house intact soil. After a plot had been 


selected, the aboveground grass (if it existed) was clipped to ground level, and 
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the litter was removed to expose bare soil before soil sampling. The PVC cylinder 
was pushed into the soil using the instruments shown in Figure 2. Firstly, the 
ground spears (8) were revolved into the ground, and the support system (7, 6) 
was set up. An air-bubble level (9) was affixed on the top of the iron roof plate 
(6) to ensure that it was level. A steel tube (1) under the roof was used to house a 
PVC cylinder (4) with a window (2) on the upper side. The internal diameter of 
the steel tube (1) was only marginally larger than the external diameter of the 
PVC cylinder to limit its shifting. A steel cap was affixed on the top of the PVC 
cylinder. Then, a vertical hydraulic jack (3) was inserted between the roof (6) and 
the steel cap through the window (2) of the steel tube (1). The PVC cylinder (4) 
was pressed into the soil by driving the hydraulic jack (3) as the roof (6) was fixed 
by the struts (7) and the ground spears (8). As the hydraulic jack (3) reached its 
maximum lifting limit, it was taken out and reset. A log (5) with the same 
diameter of the PVC cylinder (4) and a height below the maximum lifting limit 
of the hydraulic jack (3) was put on the steel cap. The hydraulic jack (3) was 
inserted into the window (2) and driven again. These operations were repeated 
until the whole PVC cylinder (4) had penetrated the soil. To reduce the resistance 
of the soil layer, the roots of plants were cut, and to avoid damage from gravel, a 
steel ring knife (10) was installed at the bottom of the PVC cylinder. The inner 
diameter of the front blade of this steel ring knife is identical to the inner diameter 
of the PVC cylinder. The inside of the upper end of the ring knife featured a 
groove, whose inner diameter equals the outer diameter of the PVC cylinder (4). 
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There were also two small vertical windows at the bottom of the steel tube (1) for 
observing whether the PVC cylinder (4) had been completely pressed into the soil. 
This instrument can smoothly and vertically push the PVC cylinder into the soil 
without shock and rotation and can thus obtain the soil core with minimal 
disturbance. After the procedure, the PVC cylinder was dragged out, covered by 
a PVC cap at each end, and wrapped in a sponge to protect the soil from 
mechanical disturbances. 
<Figure 2> 


2.2. X-ray computed tomography scanning and image analysis 


A Philips Medical Scanner (Royal Dutch Philips Electronics Ltd., Model: 
iCT 256, Amsterdam, Netherlands), a 256-slice spiral X-ray computed 
tomography (CT), was used to scan the obtained soil cores. The energy level was 
set to 140 kV and the X-ray tube current was 153 mA. The scanning was 
continuous with a slice thickness of 0.67 mm and a scan interval of 0.3 mm. After 
scanning and reorientation, an image sequence was produced with more than 
1667 slices for each soil core, in a coronal view. Voxel Cone Tracing acquired 
16-bit and 1024x1024 images with a voxel dimension of 0.122 mm x 0.122 mm 
x 0.3 mm (width x length x depth). 

Images were analyzed using ImageJ software version 1.4.3.67 to examine 
the pore characteristics according to the following steps: 

(1) Slices selection: To avoid voids and minimize disturbance caused by 


sampling, 100 slices (3 cm in depth) at the soil core head and 60 slices at the soil 
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core end were removed using the “Slice Remover" tool in ImageJ. Finally, 1500 
slices (45 cm in depth) were selected for the following analyses. 

(2) Region of interest selection: The "Region of Interest" tool was used to 
exclude voids near the core walls and to minimize the effects of beam hardening. 
A square area in the center of the soil core, with a diameter of 80 mm, was selected 
as region of interest and cropped for further analysis. 

(3) Thresholding: The attenuation value of the obtained images ranged 
from -1024 to more than 3071. The intensity value of the water phantoms (mean 
— 0) was set as threshold value to differentiate air-filled spaces from other regions 
for image analysis. Values below the threshold were identified as air-filled pores, 
while values above the threshold were identified as non-pores (Feng et al., 2003; 
Hu et al., 2015). This threshold was also used by Udawatta et al. (2008) and 
performed well in the experiments for this study. Small PVC cylinders with 
known diameters were inserted into an undisturbed soil core and scanned, and 
this threshold analysis enabled interpretation of their pore sizes. After 
thresholding, the images were translated into binary images, where black areas 
represent the soil matrix and white areas represent macropores. 

(4) Macropore analysis: After segmentation, the connectivity of soil 
macropores was estimated using the plugins “Purify” and “Connectivity” of 
ImageJ (in turn). The numbers of macropores, the macropore area, macropore 
perimeter, circularity, feret, and minimum feret of each slice were calculated 
using the "Analyze Particles" tool of ImageJ in 2D. A slice porosity was 
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calculated by the pore area of slices dividing the soil core area (5024 mm?). 
Previous studies classificated pores with equivalent pore diameter (9) > 0.1 mm 
as macropores (e.g. Harvey and Nuttle, 1995), and thus, the pores detected using 
CT (at a resolution of 0.122 mm x 0.122 mm) could all be defined as macropores. 
Other researchers reported that only pores with g > 1 mm can be divided into 
macropores (Carey et al., 2010), and therefore, the present study classificated 
pores with o > 1 mm as wide macropores to better distinguish them. 

The numbers of macropore tubes, macropore volumes, and macropore 
surfaces were calculated with the plugin *3D objects counter". These 3D pores 
were divided into wide macropores or macropores, depending on whether the 
volumes exceeded 1 mm? or not, respectively. 3D visualization of soil macropore 
networks in the soil cores was reconstructed using the plugin “3D viewer", and 
then, the skeleton of the 3D soil macropore networks was established using the 
plugin “Skeletonise 3D". The numbers and lengths of both macropore tubes and 
their branches were analyzed using the plugin *Analyse skeleton". 

«Figure 3» 


2.3. Statistical analyses 


Differences among measured parameters within-group sites were analyzed 
using one-way analysis of variance (ANOV A) combined with Fisher's protected 
least significant difference (LSD) test. All statistical analyses were conducted 
using SPSS software (version 19.0, IBM Inc., USA), applying a p = 0.05 level of 


confidence. 
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3. Results 


3.1. Soil properties 


Across the investigated catena, the soil properties varied greatly from the 
ridge to the valley (Table 1). The soils were clay loam, loam, and sandy loam at 
the ridge, slope, and valley, respectively. The sand contents increased from the 
ridge to the valley, while the clay contents followed the opposite trend, and silt 
contents were highest at the slope. The bulk densities at the ridge were 
insignificantly higher than at the slope, while bulk densities at both sites were 
significantly higher than at the valley floor for all depths. However, the 
differences of total porosity contrasted with the bulk density. The capillary 
porosities increased from the ridge to the valley but only showed significant 
differences between the ridge and the valley at 20-40 cm. No significant 
differences were found between non-capillary porosities of the 0-20 cm soil 
layers at different sites, however, at the 20-60 cm layer, the non-capillary 
porosities at the valley were significantly higher than both at the ridge and slope. 
The saturated soil moisture increased from the ridge to the valley; however, the 
field water capacity only showed a significant difference at the 20-40 cm layer. 
The aggregate weight percentages and MWD increased from the ridge to the 
valley. The K of the ridge and slope were 135.7 + 98.7 and 80.6 + 56.1 mm h”, 


respectively. Ridge sites had higher Ks than slope sites, but their differences were 
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not significant (p = 0.449). The Kps at the valley could not be measured as it 
exceeded the maximum measurement range of the DualHead Infiltrometer (1150 
mm h`’). 


3.2. Visualization of macropore networks 


Visualizations of soil macropore networks in soil columns are shown in 
Figure 4. Soils at the valley have the greatest macroporosity and more macropores 
than soils at the ridge; however, the soils at the slope have minimum 
macroporosity and macropores. At the ridge and slope, soil macropores are 
mainly distributed at the surface layer (0-15 cm), where they formed networks, 
while at deeper soil layers (>15 cm), soil macropores were fewer and isolated. 
Exceptionally, the Ridge #3 soil core showed many macropores at 15-35 cm 
depth, because this sample contained an ants nest at this depth. The soil columns 
of the valley showed high macroporosity at all depths (Valley #3). Sometimes, 
the macropores at the deeper layer (30-40 cm) showed higher numbers than those 
at the shallower layer (0-15 cm, Valley #2). The directions of soil macropores at 
the ridge and slope were mainly vertical with a few horizontal networks at the 
surface layer; however, the soil macropore directions of the valley were mainly 
horizontal, which differed from the other two sites. 


3.3. CT-measured macropores and macroporosity of whole columns 


The total number, volume, and surface of soil macropores, as well as the 
macroporosity at the valley, were significantly higher than those at the ridge and 


slope (Table 2). The wide macroporosities (volume > 1 mm?) were 3.3 + 1.9%, 
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1.4 + 1.4%, and 10.5 + 3.3% for soil columns at the ridge, slope, and valley, 
respectively, and the numbers of wide macropores (volume > 1 mm?) were 428 + 
75, 266 + 188, and 1468 + 194, respectively. Although the numbers of wide 
macropores only supplied a small part of the total macropores detected by CT, 
the large macropore volumes and wide macroporosity accounted for the majority 
of the total macropore volume and macroporosity, respectively. The mean 
volume and mean surface of macropores did not show significant differences 
between different sites. The mean connectivity at different sites showed great 
differences (Table 2); however, the differences were insignificant since there 
were great variations among each site. 


3.4. Spatial distribution of macropores 


Large variations were found in the macropore number distribution along 
the depth for different sites (Figure 5). Although the trends were different in 
different replicates, the average trend of macropore numbers decreased with 
increasing soil depth. At the top 3 cm of the soil cores, the wide macropore 
numbers were 11.4, 10.4, and 27.5 at the ridge, slope, and valley, respectively. At 
the ridge, the average wide macropore number decreased to 7.4 at a depth of 3-6 
cm and ranged between 4.8 and 3.7 at a depth of 33-45 cm. The average wide 
macropore number at the slope decreased to 3.8 at a depth of 6-9 cm and remained 
within the range between 1.9 and 3.2 at a depth of 27-45 cm. The wide macropore 
numbers at the valley were highest among the three sites with a maximum number 


of 27.7 at a depth of 3-6 cm. The numbers of total macropores also decreased 
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with increasing depth (Figure 5). The percentages of the wide macropore number 
to the total macropore number increased with soil depth at both the ridge and the 
valley. 

Soil macroporosities differed greatly among the three investigated sites 
(Table 2), and also varied along depths based on the data of slices (Figure 6). Soil 
macroporosities mainly decreased with increasing soil depth except for Ridge #3 
and Valley #2. At the top 3 cm of the analyzed soil cores, the mean 
macroporosities were 0.62%, 0.48%, and 2.04% for the ridge, slope, and valley, 
respectively. The macroporosities at the ridge remained above 0.50% over 12 cm 
depth and then became lower than their mean macroporosity of 0.34% except for 
the depths of 15-18 cm and 30-33 cm, where extended ant nests were accidentially 
sampled at Ridge #3. Ant nests are common in this region and it cannot be 
avoided that these are included in the sample. At the slope site, only 
macroporosities at soil layers above 9 cm depth exceeded 0.48%, and the 
macroporosities of the soil layers below this depth did not exceed 0.16%. The 
macroporosities in the valley were highest among the three sites for each layer, 
ranging from 0.27% to 1.33% at depths below 18 cm. The trends of the wide 
macroporosities were similar to those of the macroporosities (Figure 6) since 
wide macropores represent the majority of total macropore volume (Table 2). 


3.5. Macropore size distribution 


The macropore sizes (9) mainly ranged between 0.1-1 and 1-2 mm, 


accounting for more than 44.0%, and 32.3%, respectively (Figure 7). o > 5 mm 
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were less than 1% of the total macropores at all three sites. At the slope site, the 
percentage of 0.1-1 mm macropore was highest among the three sites, and the 
percentages of the wide macropores larger than 2 mm were lowest. The 
percentages of 1-4 mm wide macropores at the valley were highest among the 
three sites, while the percentages of wide macropores >4 mm were highest at the 


ridge. 


4. Discussion 


Macropore networks were distinctly different among different positions at 
the collapsing gullies and different types of macropores were detected by CT. 
Macropores formed by earthworms were highly continuous, relatively large, and 
tubular (Luo et al., 2010); however, earthworms were not found at the ridge and 
slope since the soil was very barren with soil organic carbon levels as low as 1.71 
g kg! (Zhang et al., 2020). However, a network of ant nests was unexpectedly 
found in the middle of Ridge #3 (Figure 4). The ant network macropores were 
characterized by several big chambers connected by tunnels, and these 
characteristics were similar to those observed before Li et al. (2017), L1 (2017). 
Macropores formed by roots were highly continuous and round in shape, the size 
of which generally decreased with depth. These macropores could be found in all 
soil cores, especially within the top 10 cm, which is in accordance with the roots 


that are mainly distributed within the topsoil layer. Smaller, more randomly 
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distributed, and less continuously macropores were likely inter-aggregate 
macropores, which formed by wetting and drying (Luo et al., 2010; Hu et al., 
2015). These are also the main macropores at soil layers deeper than 20 cm at the 
ridge and slope sites. Another phenomenon is the presence of many horizontal 
macropores at the valley site. This is because plants generally have shallow roots 
and grow horizontally (e.g., bamboo) because of the high soil water level and the 
absence of water deficits in the valley. 

In addition, the erosion and sediment process greatly affected the number, 
size, macroporosity, and continuity of soil macropores. Soil macropores at the 
slope were comparatively rare and small (Table 2), similar to previous findings 
of Zhang et al. (2016). Soil erosion was more intensive at the slope compared to 
the ridge (Sabzevari and Talebi, 2019), and the soil horizon A at the slope had 
eroded during observation. Soil horizon A usually has a higher porosity and more 
macropores than horizon B (Hu et al., 2020b). The barren slope soil without 
horizon A also restricted the growth of plants and the P. massoniana growing 
there were aged and dwarfed without undergrowth, which further limited the 
macropore development. However, at the valley, the soil was deposited by floods 
from the slope, and many dead tree branches were buried in different layers, 
which then formed macropores after decomposition. Also, the sediment mainly 
consisted of sand with an abundance of inter-aggregate macropores. 

Macropore characteristics altered at different sites, which affected the soil 
water movement and the processes of soil erosion and collapse. At the slope, few 
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macropores limited the infiltration capacity (Iversen et al., 2012; Hu et al., 2020a). 
Superimposing the flow from the upslope, the overland flow at the slope was 
strong and erosive. Severe erosion limited soil water because of few infiltrations, 
and soil barrenness caused by the loss of horizon A inhibited the establishment 
of plants. Thus, the interaction between the low number of soil macropores, low 
vegetation cover, and fierce overflow at the slope was mutually promoting. The 
development of slope erosion thus stretches into the collapsing gullies. This could 
be one of the evolutionary mechanisms of collapsing gullies after the destruction 
of the original vegetation. This agrees with previous literature, which proposed 
that collapsing gullies developed because of the loss of vegetation (Zhang and 
Zhong, 1990). 

The ridge had more macropores, which were larger in size and volume, and 
the macroporosity was also higher than that of the slope. This was caused by the 
relatively weaker soil erosion and better vegetation conditions at the ridge 
compared with those at the slope. Activities of plant roots and soil macrofauna 
(e.g., ants) are advantageous for the formation of macropores. These conditions 
promote water infiltration at the ridge. According to former research, the soil 
water content is an important factor that promotes the development of an existing 
collapsing gully (Shi, 1984; Li, 1992). More rainfall infiltration decreases the 
shear resistance of the soil layer while it also adds to the weight of the soil column. 
These developments further increase the probability and risk of soil collapse and 
the development of a collapsing gully. Macropores decreased with increasing soil 
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depth, indicating that the infiltration ability decreases with increasing soil depth 
as Duan et al. (2018) reported. Thus, the soil water would remain at the topsoil 
layer, thus increasing the weight of the soil layer and decreasing the shear 
resistance of soil (Deng et al., 2018). 

The valley site was covered by sediment from the upslope and had an 
abundance of macropores. Consequently, infiltration would be large and the 
underground flow may be strong under these conditions. The high soil water 
content and subsurface flow increase the risk of subsurface flow erosion (Tebebu 
et al., 2010). The sediment was sandy, containing only a little clay. It has been 
reported that alluvial sediment is very erosive when the vegetation cover is 
insufficient (Jiang et al., 2020). As observed in the field, once the check dam was 
completely filled, the overland flow could easily erode the inside sediment and 
thus form a deep gully. Once the sediment was removed, the base of the 
collapsing wall would be lost, and the risk of collapse would increase without 


base support. 


5. Conclusions 


The purpose of this study was to quantify the macropore features of soil 
using CT at different sites of a region with collapsing gullies. The results yielded 
macropore information about the effects of soil erosion and sediments on the soil 


structure and clarified an aspect of their interaction. The ridge had more and larger 
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macropores than the slope, while the sediment at the valley had the most 
macropores. The numbers of macropore and their macroporosities mostly 
decreased with increasing soil depth. The equivalent pore diameter of macropores 
mainly remained below 2 mm. These characteristics differed according to the 
tomography position, soil erosion, plant cover, and soil macrofauna. The high 
macroporosities of the valley and ridge increase the risks of subsurface flow 
erosion and soil collapse, which would promote the development of collapsing 
gullies. The soil macropores at the slope were rare and infiltration abilities were 
comparatively small; thus, overland flow and slope erosion can be assumed to be 


strengthened. These effects contribute to the formation of collapsing gullies. 
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Table 1 Soil properties at the ridge, slope, and valley of collapsing gully region. Values are average + standard deviation, based on three replicates. 


Capillary Non-capillary Saturated soil Field water Aggregate weight 


Sites Depth Sand Silt Clay Bulk density Total porosity MWD 
porosity porosity moisture capacity percentage 
cm gkg! gkg! gkg! g/cm? 96 % % % (W/DW) 926(W/DW) %(W/DW) mm 
0-20 367.74£75.5a 312.3419.8 320.0457.4a 1.49+0.04a 40.48x1.48a 37.99+2.56a 1.65+0.06a 27.2441.53a 21.11+1.98 28.13+7.44a 1.03+0.25a 
Ridge 20-40 1.47+0.04a 41.9622.19a 39.6043.40a 1.40+0.25a 28.67+42.36a 22.7843.15a 31.13414.74a 0.87+0.2a 
40-60 1.47+0.03a 41.8222.82a 39.2744.14a 1.49+0.24a 28.4542.12a 22.7543.50 25.53410.24a 0.7430.1a 
0-20 464.6x16.0ab  409.4+83.7 126.0£100.6b — 1.3150.11a 43.1742.20a 41.3942.2ab 1.79+0.28aA 33.0742.44a 25.02+4.57 34.99+17.89a 1.12+0.18a 
Slope 20-40 1.41+0.04a 43.66+2.06a 42.1122.03a 1.54+0.09aAB 30.9842.19a 24.54+3.84a 37.82+£12.25b 0.98+0.14b 
40-60 1.42+0.02a 43.38x1.41ab 41.96x1.44ab 1.42+0.06aB 30.65+1.42a 24.91+3.02 40.45+8.85b 1.01+0.1b 
0-20 545.3+60.4b 335.8+44.8 118.9+15.7b 1.08+0.14b 49.46+6.13b 47.50+6.46b 1.95+0.36a 46.65+10.87b 29.94+10.3 65.15+27.63b 1.65+0.37b 
Valley — 20-40 0.98+0.04b 51.55+0.94b 49.34+0.51b 2.22+0.61b 52.46+2.93b 32.76x3.92b 45.7X11.79bc 1.3120.32c 
40-60 1.1120.13b 46.7243.03b 44.7542.95b 1.97+0.18b 42.50+7.43b 21.74+2.68 47.73x17.95b 1.53+0.29c 


MWD: mean weight diameter of aggregates; W: weight; DW: dry weight. At the same depths, sites with different lowercase 
letters indicate significant differences (p < 0.05), while the same lowercase letter indicate no significant differences (p > 0.05). 
At the same sites, Layers of different depth with different capital letters indicate significant differences (p < 0.05), while the 


same capital letter indicate significant differences (p > 0.05). 
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617 Table 2 Characteristics of macropores of different stages of the slope. V is the volume of macropores. 


Stages Macropore number Volume Surface Macroporosity Mean macropore size Mean surface Connectivity 
Total V2imm? Total V2imm? Total V>1mm’ Total V>imm V2>1 mm V> 1mm? 
cm? cm? m? m? % % mm? mm? 
Ridge mid-erosion 1828+324b 428+75b 7.8X4.4b  7.5+4.3b 2.941.2b  2.6+1.1b 0.34+0.19b 0.3340.19b 16.8+7.4 58.8x18.5 15.1+25.6 
Slope max-erosion 1189+747b 2664188b 3.443.2b 3.2+43.1b 1.741.4b 1.541.3b 0.15x0.14b  0.14+0.14b 10.6+2.9 51.8+8.6 5.948.9 
Valley deposit 5919+703a 1468+194a 24.747.5a 23.747.5a 10.4+42.6a 9.2+2.6a 1.0940.33a 1.0540.33a  15.943.2 62.1+9.8 42.3+40.1 
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Figure cations: 

Figure 1 (a) The location of soil core sampling sites in Guangdong Province, China, and 
(b) their geomorphic and environment around. Shaded provinces in (a) represent 
the main distribution areas of Benggang. 


Figure 2 The instrument used to push the polyvinyl chloride cylinders into ground. 1 - steel 
tube, 2 - window, 3 - lifting assembly (vertical hydraulic jack); 4 - sampling tube (PVC 
cylinder); 5 - log plate; 6 - steel roof; 7 - struts; 8 - ground spear; 81 - ring; 82 - limit part; 83 - 
crack; 84 - down warping part; 85 - up warping part; 9 - level; 10 - steel ring knife; 101 - front 
blade; 102 - groove. 

Figure 3 Procedures of image analysis in the study. 

Figure 4 Three-dimensional visualization of soil macropore networks in the soil columns and 
soil macroporosity along the soil depth at the ridge, slope, and valley. The number after ridge, 
slope, or valley is the replicated column number. 

Figure 5 The macropore numbers along the soil depth at (a) the ridge, (b) the slope, and (c) the 
valley. The black, red, and green solid lines are the macropore numbers of three replicates with 
the blue solid line being their mean for each 3 cm depth. The blue short dash lines are the 
number of macropores with equivalent pore diameter » 1 mm. 

Figure 6 Distribution of macroporosity along the soil column depth at the ridge, slope, and 
valley. 


Figure 7 The frequency of equivalent pore diameters of soil columns at different sites. 


29 


chinaXiv 


Ridge site SM ESE cS 


Valley-site 


30° 


10° 


* Soil cores sampling site 
0 500 1,000 2,000 km 
L— o nnd 


Figure 1 (a) The location of soil core sampling sites in Guangdong Province, China, and (b) their geomorphic and 


environment around. Shaded provinces in (a) represent the main distribution areas of Benggang. 
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Figure 2 The instrument used to push the polyvinyl chloride cylinders into ground. 1 - steel tube, 2 - window, 3 - 
lifting assembly (vertical hydraulic jack); 4 - sampling tube (PVC cylinder); 5 - log plate; 6 - steel roof; 7 - struts; 
8 - ground spear; 81 - ring; 82 - limit part; 83 - crack; 84 - down warping part; 85 - up warping part; 9 - level; 10 


- steel ring knife; 101 - front blade; 102 - groove. 


31 


650 


651 
652 


Image 
acquisition 


Slices selection 
and ROIs crop 


Image 
thresholding 


v 


Macropore 
analysis 


Slice 


Macropore numbers || Macropore area 


Macroporosity | Macropore perimeter 


Circularity 


Feret 


Minimum feret 


v 


Column 


Numbers of macropore tubes 


Macropore volumes 


Macropore surfaces 


Connectivity 


Rendering and 
visualization 


i 


Macropore 
analysis 


>| Column 
L 


Numbers of macropore tubes and their branches 


Lengths of macropore tubes and their branches 


Figure 3 Procedures of image analysis in the study. 
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= 656 Figure 4 Three-dimensional visualization of soil macropore networks in the soil columns and soil macroporosity 
along the soil depth at the ridge, slope, and valley. The number after ridge, slope, or valley is the replicated column 


number. 
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Figure 5 The macropore numbers along the soil depth at (a) the ridge, (b) the slope, and (c) the valley. The black, 


red, and green solid lines are the macropore numbers of three replicates with the blue solid line being their mean 


for each 3 cm depth. The blue short dash lines are the number of macropores with equivalent pore diameter > 1 


mm. 
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Figure 6 Distribution of macroporosity along the soil column depth at the ridge, slope, and valley. The 9 is 


668 equivalent pore diameter. 
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671 Figure 7 The frequency of equivalent pore diameters of soil columns at different sites. 
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